AC- and DC-driven noise and I-V characteristics of magnetic nanostructures 
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We study a structure consisting of a ferromagnetic (F) layer coupled to two normal metal (N) 
leads. The system is driven out of equilibrium by the simultaneous application of external dc and ac 
voltages across the N/F/N structure. Using the Keldysh diagrammatic approach, and modeling the 
ferromagnet as a classical spin of size S 3> 1, we derive the Langevin equation for the magnetization 
dynamics and calculate the noise correlator. We find that the noise has an explicit frequency depen- 
dence in addition to depending on the characteristics of the ac and dc drive. Further, we calculate 
the current-voltage characteristics of the structure to O [l/S^) and find that the nonequilibrium 
dynamics of the ferromagnetic layer gives rise to corrections to the current that are both linear and 
nonlinear in voltage. 

PACS numbers: 72.25.-b; 75.70.Cn; 75.75.-c 



I. INTRODUCTION 

Magnetization dynamics in small nanomagnets 
has recently attracted a lot of theoretical^"— and 
experimental^'^ attention due to advances in man- 
ufacturing magnetic nanostructures. The topic of 
magnetization noise has become an exciting subject ow- 
ing to its possible influence on magnetization switching^ 
and conductivity of these structures 4 It has been shown 
that the noise in magnetic structures, such as spin valves, 
can be colored, i.e., it can have a nontrivial frequency 
dependencei^ In diffusive metallic conductors colored 
noise has been observed experimentally,^- however in 
magnetic structures it still requires further investigation. 

Experiments involving magnetic nanostructures typi- 
cally involve the simultaneous application of dc and ac 
voltage where the ac bias is found to aid the magnetiza- 
tion switching. Therefore in this paper we study a nor- 
mal metal/ferromagnet /normal metal (N/F/N) structure 
which has been driven out of equilibrium by the simul- 
taneous application of a dc and ac voltage. We show 
that the effect of this driving is to produce a noise in the 
magnetization dynamics that is colored. In addition we 
determine how the I-V characteristics of the device are 
affected by the dynamics of the ferromagnetic layer. 

The schematic of the N/F/N structure we study is 
shown in Fig. [TJ The ferromagnetic layer is assumed to 
be small so that it may be modeled as a single-domain 
magnet. At the same time, the spin of the magnet is con- 
sidered to be large (spin S 1) so that it can be treated 
as a classical variable. The aim of this paper is twofold, 
one is to derive the Langevin equation for the magneti- 
zation dynamics, and second is to present a calculation 
of the I-V characteristics. In the absence of any magne- 
tization dynamics, the N/F/N structure is Ohmic,^^ We 
show that the dynamics of the ferromagnet gives rise to 
corrections to the I-V characteristics that are both linear 
and nonlinear in voltage. 

The paper is organized as follows. In section |TT] we 
present the model. In section Hill we study the nonequi- 
librium properties of the model in the limit 5 — >■ oo. 



when the magnetization is static. In section |TV] we study 
small fluctuations of the magnetization about the order- 
ing direction thus deriving the Langevin equation and the 
noise spectrum. The results of this section are then used 
in section IVl to calculate the corrections to the current- 
voltage characteristics arising due to the magnetization 
fluctuations. Finally in section [VI] we summarize our re- 
sults. 



II. MODEL 

We consider a model Hamiltonian H = iJ.,„ + IIi+ Ht, 
where Hm describes the ferromagnetic layer, Hi repre- 
sents the two normal-metal leads, and Ht models the 
tunneling between the leads and magnetic layer. The 
Hamiltonian for the magnetic layer Hm is 



= -{DSl + BS,) + J ^ S • s, 



kCT 



[dij^... (1) 



Here the first term models a material (or shape) 
anisotropy with the anisotropy constant the second 
term describes the interaction of the macrospin S with 
the magnetic field B applied for simplicity in the same z- 
direction as the anisotropy. The third term describes the 
interaction of the macrospin with the spins of the itiner- 
ant electrons as in the s-d model. It can be rewritten as 
JY.i{SzSzi + S+S-i-^ S-S+i), where S± = {Sx±iSy)/2, 



S±i 



and 



with Gaji being the components of Pauli matrices. Here 

(ijj.^(c?kcr) creates (destroys) an electron in state with mo- 
mentum k and spin component a. The pure macrospin 
part of the Hamiltonian can be rewritten as —{DS^ -\- 
BSz) — const -)- bS-i-S-, where the constant part is 
- BS and b = 4(0 + B/{2S)). Nanomagnets are 
typically characterized by a significant anisotropy. This 
along with the fact that S ^ 1 implies that the fluctua- 
tions of the nanomagnet about the ordering direction are 
small. Our theoretical treatment will therefore involve a 
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FIG. 1: (Color online) A sketch of the N/F/N structure. The 
left and right normal-metal leads are coupled through tunnel 
barriers to the ferromagnetic layer. 



perturbative expansion in spin fluctuations, which as we 
shall show is equivalent to an expansion in l/S. 

We assume that the electrons in the leads are non- 
interacting. To model the ac bias voltage we introduce 
a time-dependence of the lead single-particle energies^^ 
namely, eka{t) = CfeQ + K2cCOs(wo< + fa) where a labels 
the left [L) or right [R) lead, and Vac, "-^o are, respectively, 
the amplitude and frequency of the ac bias. Thus, the 
lead Hamiltonian is 



V 

Cko 



(2) 



ka.a^L.R 



The coupling between the leads and magnetic layer is 

Ht^ E (^"4,„.„dk.+H.c.). (3) 

]i.,ko.a,a^L,R 

In Eq. ^ k ~ (k, ko) where k is a two-dimensional mo- 
mentum in the plane perpendicular to the tunneling di- 
rection and is assumed to be conserved on tunneling. 

To study this nonequilibrium problem we employ the 
Keldysh formalism.'^'-'^'^ We introduce variables S"^' = 
(S+-|-S-)/2and S« = (S" -S+)/2 where the upper ± in- 
dices correspond to the time-ordered (anti-time-ordered) 
directions on the Keldysh contour, and the Keldysh path 
integral takes the form, Zk = / 2?[S'=', S«] e-"^^' . Here 
Sk is the effective action for the macrospin obtained for- 
mally by integrating out all fermionic degrees of freedom: 



Sk ^ 26Tr {S^^Sl + SlS^_') 



+ «Tr In 



5da 



} — — JSz—^ — JS-L<7_ — JS^a, 



,(4) 



where Sa^ 



■ z,± 



-S and gda 



Qq Qcl ] ^^^^ 'jaa l r. A 
J V ^ yd<7 

the Green's functions of the free electrons in the mag- 
netic layer, cr^ and a± = {a^ ± (Ty)/2 are Pauli matrices, 

and ^ — ( Q j^A ] is the self-energy due to coupling 

to the leads. As we shall show, E depends on the ac and 
dc bias and is independent of a because the leads are 
non-magnetic. In what follows we make the assumption 
that the fluctuations of the macrospin from the ordering 
direction are small. Thus we write JSz = JS + J5Sz, 
and eventually expand Eq. ^ perturbatively in the fluc- 
tuations 5Sz , S± . 



III. GREEN'S FUNCTIONS IN THE 
MAGNETIC LAYER 

We first discuss the properties of the nonequilibrium 
system when the magnetization does not fluctuate. De- 
noting Go to be the Green's function of the electrons in 
the magnetic layer when 5Sz = S± =Q, Eq. (g]) implies 

[GlY' = K]-^-S^-J5|, (5a) 

(5b) 



The symbol o in Eq. ([5]) denotes convolution in the time 
domain, and the self-energies due to coupling to leads are 



I]«(-)(M') = E^^.9f^"'(M'). 

feo.Q 



(6) 



dkaa'^ are the retarded and Keldysh components of the 
electron Green's function in the leads and are defined as 

9ka<.it,t') = -*([cfe.a(t),cL„(t')]), (7b) 



e-'/!»<**i^'==(*i)cfc^„(-oo), we find 



Since Ckaa{t) 



(8a) 



gtait,t') = - 2/(.L - /^.)]e-''^*' (8b) 

where / is the Fermi distribution function in the leads 
and we have used that ([cfeo-Q (— oo), cj^^^(— oo)]) = 1 — 
2/(^fcQ "Ma) where is the chemical potential of lead 
a and a dc bias corresponds to /iL 7^ fiR. Using the 

identity e*(°"") = Z]J^-oo ""•^n(^)' where Jn{z) are 
Bessel functions of the first kind, we find 

5^_(i,t') = -*^(^-t')e-<-'*-*'^ 



Vac 

UJQ 



intjQ {t — t' ) 



——00 



^iifia (m— n) —iujQ irnt' — nt) 



(9) 



Changing variables to t — t — t' and T — {t + 1')/2 one 
may write mf — nt = (n — m)T + (n -\- m)T /2. In what 
follows we average Green's functions over time T ^ i^t^^ , 
where the averaging is denoted by gka{t,t'). This is 
justified when the magnetization dynamics is slow com- 
pared to ujq^ (a precise condition for this will be given 
in SecHVl). 

As a result of the time averaging, terms corresponding 
to rt 7^ m in Eq. ^ vanish. This leads to a time-averaged 
retarded self-energy S-^(a;) = ^*X]aeL B.^a where Fq = 
TTi't'^ is the decay rate of the electrons into the leads, and 
1/ is the density of states in the leads. In what follows 
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FIG. 2: (Color online) Distribution function of electrons in the 
magnetic layer /„eq(aj) given by Eq. (a) pure ac voltage 

case, fj.L ~ fJ-R ~ and Vkc/i^o = 0.5; (b) the case of nonzero 
ac and dc voltages, the parameters are Vl = Tij, fiL = — A'fl ~ 
1.5a;o, and 14c /i^o = 0.5. Frequency uj is measured in units of 

LUQ. 



we will assume Fq to be independent of energy. From 
Eq. ([6]) , the time-averaged Keldysh component of the self- 
energy becomes 



-2i^F„ 



n— — OO 



(10) 

The above discussion implies that the spectral func- 
tion of electrons in the magnetic layer Aqct = — Ini[G^o.] 
is the same as in equilibrium, Aocr(k, w) = r/[(a; — e^ — 
ctA)^ -I- F^], where F = F^ -I- F^ and the exchange split- 
ting A = JS/2. We win assume A > F so that the 
ferromagnetism of the conduction electrons is well de- 
fined. Further, the nonequilibrium distribution function 
fneq of the clcctrons in the magnetic layer (defined as 
2iAo^(k,a;)[l - 2/„eg(t^)]) is 



Vac 

OJQ 



/(w - nwo - AIq), (11) 



Below we consider the case of zero temperature when 
the Fermi function /(w) = 0{—uj). A typical fneqi^) is 
plotted in Fig. [21 While in the pure dc case fneqiy^) is 
a weighted sum of Fermi functions of the left and right 
leads, the ac bias adds steps to fneq at frequencies ±|n|a;o 
corresponding to photon absorption and emission. 



IV. LANGEVIN EQUATION 

We expand the effective action (U) to quadratic order 
in the fluctuations to obtain, 

Sk = 2bTr{SfS'L + SlS^_') + {SfU^.S^ + H.c.) 
1 



where 11^^ are the components of the polarization oper- 
ator that are calculated following standard techniques.— 
Note that Ilxxit^) — ^yyi^^) and Ilxyiuj) = —Ilyx{uj). 
Moreover, to leading order in spin fluctuations II^^ does 
not play a role. For small frequencies w <C A, we find: 



^xxi^) - -^l^xxUJ, n^'Juj) ~ -i(3xyUJ, 



/3x 



Jxxl^, 

F 



A2 + F2 ' 



^xy 



xyy 

A 

A2 + r2 ■ 



(13a) 
(13b) 



Note that /S^y = yPxx- 

The action, Eq. p^ . may be diagonalized in the basis 
5^' — {S^±iSy)/2 thus yielding the Langevin equation^ 



{l3xxTilixy)Si^U. (14) 



where — (^^ ± i£^y)/2 is an auxiliary field 

representing noise whose correlator is given by 
{^a=x,y{^)^b=x,y{-i^)) = «n^(w). Wc havc found an 
analytical expression for 11^ in terms of a double sum 
over squares of Bessel functions. For Vac ^ we may 
keep only terms corresponding to single photon absorp- 
tion and emission processes. For = —^J■R — V/2 and 
a; ^ A the noise correlator is 



{Cx{^)^xi~^^)) 



Wx 



F2 



^^xxi^) 

{\uj + V\ + \lj-V\ 
F- 



r2 

F2 



I 2wo 



E 



F2 



(15) 



Similarly, the off-diagonal component (f2;(w)^j,(— w)) = 
i^xyi^) is 

AFCJ ^r- .9 fVac\ j2 fVac 

Wo 



(A2 +F2) 
EqF^ 



Wo 



^ r2 

a,l3=L,R 



Icj 4- (m - n)uJo + fJ^a- fJ^p] -(16) 



An effective temperature T^e may be extracted from 
the zero frequency limit of the noise correlator. 
Eq. CS]) implies Teff - (TLTR/r^)V + (Kc/2wo)2[wo + 
{2TlTr/T^){V - ujQ)e{V - Wo)] and therefore has a dis- 
continuityi^ at V = wo. In the opposite limit of w 3> A 
the noise correlator vanishes as ~ 1/w as expected. 

Equation (fT4|) can be rewritten in the form of a stochas- 
tic Landau-Lifshitz-Gilbert equatiorti^, 



S = 7S X Hoff - aoS X S + 



(17) 



where 7 is the gyromagnetic ratio, Hcff = bz/{'-ff3xy) is 
an effective magnetic field, the noise is f' — tt— z x f , and 

Pxy 

the Gilbert damping constant is ao = Pxxl Pxy = E/A. 

In order to determine how the magnetization dynam- 
ics affects the I-V characteristics we will need the spin 
response and correlation functions. From Eq. ()12p the 
spin-spin response function is 



D\{u) = ^{S^l{u)Sl{-u)) = 



1 



b - PxyiO - iPxxUJ 

whereas the spin-spin correlation function is 



(18) 



i?^+(w)=»(5!_'(w)5f(-w)) 



(n^;-*nfj(w) 



\b - PxyUJ 



il3xx^\ 
(19) 
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The exchange contribution to the self-energy is 



FIG. 3: (Color online) The diagrams for one-loop corrections 
to the electron Green's function in the magnetic layer due to 
coupling to spin fluctuations: (a) exchange (Fock) contribu- 
tion E^^, (b) Hartree contribution Sjjf . Wavy hnes corre- 
spond to spin-spin correlators (5_(aj)S+(— oj)). 



and D^''_^\A) ~ Z)"'^''(— A). As expected in equilib- 
rium (1/ — Vac — 0) the components of both D(uj) and 
n(w) satisfy the fluctuation-dissipation theorem. 

It is instructive to take the inverse Fourier transform of 
Eq. (IT^ to obtain the time dependence of the transverse 
spin-spin correlation function, 



{S^_^{t > 0)5f (0)) = 



e ^ ^1 



(20) 



where wi = b{p^y-ip^^) / {Pl^ pl^), t ^ vJ^/{bT), and 
Ti = 2vJ/{hS). Equation (PO)) shows that the spin-spin 
correlations decay with the characteristic time r. Thus 
as long as I/wq ^ t, the macrospin dynamics is rather 
slow and we can use the time-averaging procedure for the 
Green's functions outlined in Sec. IIIIl 



iJ"^ f dtu 



I]|«(k,17) = _ / [G«(k,a. + J7)i^5+H 



2 J 2tt^ 



+G^J\<L,Lo + n)D^+{uj) 



(23) 



Keeping terms to leading order in J^, S*^^ is purely real 
and given by 



J' 



s|^(k,^^) = -^E 



Ea 

E 



— + arctan ^ ^ — ^ 



f Vac 

V2wo 



7r+ ^ 



arctan 



+ A + mcjQ — ^J'a 



m=±l 



(24) 



The Hartree contribution to the self-energy is given by 

Ejf(k,n) = -^Dt(--0) j vde' j ^G^,ie',co') 

where D^_{uj = 0) = 1/6 and we have set ^j., / i/de' . 
Note that the Hartree contribution is independent of 
external frequency and momentum, and therefore only 
shifts the position of the pole of but does not con- 
tribute to the corrections to the current. 

We denote the total current averaged over time T ^ 
1 /wo as I = Iq + 61 where 



n 1 



(25) 



V. CURRENT- VOLTAGE CHARACTERISTICS 

We will now study how the current- voltage character- 
istics of the magnetic junction are affected by the mag- 
netization dynamics of the ferromagnetic layer. We em- 
ploy the Jauho-Meir-Wingreen formula^^ for the tunnel- 
ing current 



k,(T 



(21) 

In the following we calculate the leading correction to the 
spectral function = — Im[G^^] due to coupling to spin 
fluctuations. The spectral function is determined from 

the Dyson equation [G^'^] = [G^q] — S^^;, where 
'Sa-d is the self-energy due to coupling to spin fluctuations. 
To one-loop order = Ej^^-f S^f , where T,%^ and E^f 
are, respectively, the exchange and Hartree contributions 
to the self-energy, see Eig. |3l To leading order in the 
fluctuations, it suffices to do perturbation theory in 



so that Gf^ 



5Gf, with 



5Gt 



riR s^R nR 



(22) 



Note that Y^^d and Ej^^ are related by A o —A. 



is the current for a static ferromagnet while 51 is the lead- 
ing correction due to spin fluctuations computed from 
Eqs. (I2I1), (1211), and dMl) for Ml = "Mfl = V/2 



SI _ TJ^ 
lo ~ 27r6A2 (1 



p2 • 



y2 

* ac 



3A2 



1/2 

ac 

2^0 



(%/5A2 



(3A2-E2)1/2 



12(A2 + E2)2 

r2)2]/2 



(A2 -hE2)2 



2 (A2 -t- E2)' 



(26) 



Since A ^ JS, this correction to the current is 0{T/bS^ ). 
Thus our perturbative treatment in spin fluctuations is 
valid as long as 6 7^ and S ^ 1. Moreover the correction 
SI > 0. This is because scattering off spin fluctuations in 
this geometry produces additional channels for electron 
conduction (in contrast to a bulk geometry where this 
scattering would cause the conductivity to decrease). It is 
worth mentioning that ac bias contributes to the Ohmic 
corrections as well with terms such as ~ (V^^/ A'^)V and 
{V^^/ujq)V. Eor the pure dc case (Vac = 0), we flnd that 
the differential conductance g = edI/dV for A ^ E is 



= 4 



e2 i/EiEfl 



1 



27r6A2 



1 



3V' 
4A2 



(27) 



Note that the quadratic in voltage corrections in Eq. 
(1271) are similar in spirit to temperature corrections (~ 



5 



T^) to the conductance. Also, our result is for a particu- 
lar choice of chemical potentials fiL = — A*-R = V/2. The 
answer in general will change if a different choice, such as 
= V and = 0, were used. The reason for this dif- 
ference is that the symmetric combination of chemical po- 
tentials (pL + (J-r) /2 plays the role of a mean chemical po- 
tential for the electrons in the nanomagnet, tuning which 
modifies the equilibrium spectral density and hence the 
linear-response conductance as well as other equilibrium 
properties. In an experiment this mean chemical poten- 
tial may be tuned by an external gate voltage. Our purely 
antisymmetric combination /i^ = ^^J■R avoids these in- 
trinsically equilibrium effects. 



tion dynamics is characterized by a frequency dependent 
noise, Eq. (fT5|) . We have also computed corrections to the 
I-V characteristics to leading (l/S"^) order in the spin- 
fluctuations. These fluctuations are found to not only 
modify the Ohmic part of the I-V characteristics, but to 
also give rise to corrections that are non-linear in volt- 
age, Eq. ([26|) . Experiments often exhibit non- linear I-V 
curves,® but the origin of the nonlinearities is usually not 
clear. The usefulness of our result is that the current is 
a function of three independent experimentally tunable 
parameters (the dc bias V , the ac amplitude Vac and fre- 
quency Wo) which can in principle allow one to extract 
the physics arising only from magnetization dynamics. 



VI. SUMMARY 
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